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Abstract. Bernstein processes are Brownian diffusions that appear in Eu- 
clidean Quantum Mechanics. The consideration of the symmetries of the as- 
sociated Hamilton-Jacobi-Bcllman equation allows one to obtain various re- 
lations between stochastic processes (Lescot-Zambrini, Progress in Proba- 
bility, vols 58 and 59). More recently it has appeared that each one-factor 
affine interest rate model (in the sense of Leblanc-Scaillet) could be described 
using such a Bernstein process. 
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1. Introduction 

The relationship between Bernstein processes and Mathematical Finance was 
first glimpsed in §5 of [9], where the Alili-Patie family of transformations S a ^ 
(see PP,§2, for the particular case a = 1, and 9., pp. 60-61, for the general case) 
was reinterpreted in the context of isovectors for the (backward) heat equation. 
Furthermore, Patie's composition formula ([S], p. 60) for these transformations was 
derived from the commutation relations of the canonical basis of the aforementioned 
isovector algebra. 

As the motivation for Patie's work had been to compute certain option prices 
in the framework of an affine interest rate model (see [9], pp. 101-104), it was 
natural to look for a direct parametrization of such a model by a Bernstein process. 
It turns out that, for any one-factor interest rate model in the sense of Leblanc 
and Scaillet ([I], p. 351), the associated square root process coincides, up to some 
(possibly infinite) random time, with a Bernstein process (§3, Theorem 3.4). The 

potential — ^ (that had been considered in ED, p. 220, on physical grounds) appears 
T 

here naturally (Theorem 3.4) : this is one more example of the deep links between 
Euclidean Quantum Mechanics and Mathematical Finance. 

2. An isovector calculation 

We shall work within the context of jS], §3. We consider the Hamilton-Jacobi- 
Bellman equation 

„, A ds e 2 d 2 s i,as\ a 

with potential 

V(t,q) = % + Dq 2 , 
q 

where we have replaced (as explained in [5]) y/h by 9. 



Date: November 17th, 2010. 



1 



2 



PAUL LESCOT 



In order to determine the Lie algebra 7-Lv of pure isovectors for (T-LJB V ), i.e. 
the algebra H of |8], we have to solve the auxiliary equation (3.29) from [8], p. 215, 
i.e. : 



, 



-\f m 2 -^q + a+\T N q{-^ + 2Dq) + l{-^ + 2Dq)+T N & + Dq 2 )- e —T N 
4 2 q 6 q 6 q z 4 

that is : 

q 2 (2DT N - - f N ) + q(-l + 2DZ) + o - Q —f N - ^ = . 
4 4 

As Tjv, Z and a depend only upon t, the system is equivalent to : 

2C7 = 

I = 2DI 

. 9 2 •■ 
o = -T N 

[ f N = 8DT N . 

Two different cases now appear : 
1) C + 

Then one must have I = 0, in which case the second condition holds automati- 
cally, and the system reduces itself to : 



a=-T N 
f N = 8DT N 



l)a) D > 

Setting e = y/8D, we find 



whence 



and 



T N = C ie ct + C 2 e~ 



J-n — — e e + O3 



therefore 



a = —T N + C4 



4 4 

where (Cj)i<j<4 denote arbitrary (real) constants. In particular 

dim(U v ) = 4 . 

l)b)D = 

Then from Tn = follows 

T N = dt 2 + C 2 t + C 3 
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for constants C\ , C 2 , C3 . Then 



2 •• 9 2 d 
a = —T N = 



4 2 
and 

Therefore, here too, dim(Hv) = 4; furthermore, we get an explicit expression 
for the isovectors : 

N t =T N = dt 2 + C 2 t + C 3 , 

N q = \qT N + l 

= ^(2dt + C 2 ) 



2 



2 



and 



7V S = -c 



1 o •■ 

= -jQ T N - ql + a 

1 2 2 
= --<? 2 .2Ci + yCi< + C 4 

A canonical basis for Wy is thus given by (Mi)i<i<4, where Mj is characterized 
by Cj — Sij (Kronecker's symbol). Using the notation of [7J, it appears that 

Mi = ^N 6 , 

M 2 = ^N 4 , 
M 3 = JVi , 

and 

m 4 = -±n 3 , 

therefore Hv is generated by iVi, N 3 , N 4 and A^. We thereby recover the result 
of [5], p. 220, modulo the correction of a misprint. This list ties in nicely with the 
symmetry properties of certain diffusions related to Bessel processes (see [6] for a 
detailed explanation). 

l)c)D<0 

Setting now e = \/—8D, we find 



whence 



T N = Ci cos (et) + C 2 sin (et) , 
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and, as above : 



C C 
Tjv = — sin (et) cos (et) + C3 , 



therefore 



6 2 C X 



cos (et) 



2 G 



2 sin (et) + C\ 



4 v ' 4 
where (C , j)i< :) <4 denote arbitrary (real) constants. In particular 

dim(H v ) = 4 . 

2)C = 

Then the system becomes 



I = 2DI 

e 2 ■• 
ff = T TAr 



T w - 8DT, 



N 



The equation for I on the one hand, and the system for (a, Tjv) on the other 
hand, are independent, and, as above, the first one has a two-dimensional space of 
solutions and the second one a four-dimensional space of solutions, i.e. 



dim(%v) = 6 . 



Whence 



Theorem 2.1. The isovector algebra Hv associated with V has dimension 6 if and 
only if C = ; in the opposite case, it has dimension 4. 



ON AFFINE INTEREST RATE MODELS 



5 



3. PARAMETRIZATION OF A ONE-FACTOR AFFINE MODEL 

As general references we shall use, concerning Bernstein processes, our recent 
survey ([S]), and, concerning affine models, Henon's PhD thesis([3]) as well as 
Leblanc and Scaillet's seminal paper ([!]). 

An one-factor affine interest rate model is characterized by the instantaneous 
rate r(t), satisfying the following stochastic differential equation : 



(3.1) dr(t) = y/ar(t) + /3 dw{t) +(<j>- Xr(t)) dt 

under the risk-neutral probability Q (a = corresponds to the so-called Vasicek 
model, and (3 — corresponds to the Cox-Ingersoll-Ross model ; cf. [J]). 
Assuming a > 0, let us set 

<P =def 9 H , 

a 

—def , 

a 

and let us also assume that 4> > 0. 

The following two quantities will play an important role : 



and 



Let us set X t = ar(t) + /9 '. 
Proposition 3.1. Let rp € R ; then the stochastic differential equation 



(3.2) dr{t) = y/\ar{t) + J\ dw(t) + (4> - Xr(t)) dt 

has a unique strong solution such that r(0) = ro- Furthermore, in case that 

ar + (3 > , 

one has ar(t) + /3 > for all t > ; in particular, r(t) satisfies (3.1). 

Proof. Let us set X t = ar{i) + f3 ; it is easy to see that, in terms of X t , equation 
3.2 becomes: 



dX t = adr(t) 



= a(s/\Xt\dw(t) + {(/>- A— — -)dt) 



= ay/\Xt\dw{t) + (a4> - XX t )dt . 

We are therefore in the situation of (1), p. 313, in [2], with c = a, a = a<j>, and 
b = — A ; the result follows. 

In case A 7^ 0, one may also refer to [3],p.55, Proposition 12.1, with a — a, k = A 

and a = ——. U 
A 
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We shall henceforth assume all of the hypotheses of Proposition 3.1 to be satis- 
fied. 



Corollary 3.2. One has 



f X t = e- At Y( a (e 4A ' 1) ) for A ? , and 
\ X t = Y{*g) for A = 

where Y is a BESQ 6 (squared Bessel process with parameter 5 ) having initial value 
Y = ar + P. 

Proof. In case A ^ 0, one applies the result of [3], p. 314. For A = 0, let 

Zt '■= —^X t ; 

cr 

it appears that 

dZ t = 2y/\Z^\dw{t) + 5dt , 
whence Z t is a S.BS'Q^process. As 

X t := *2k , 
t 4 r ' 

the scaling property of Bessel processes yields the result. □ 
Theorem 3.3. If 6 > 2 one has, almost surely : 

Vt > X t > ; 

on the other hand, if S < 2, almost surely there is at > such that X t = 0. 
Proof. We apply Corollary 1, p. 317, from [4] (12.2) yielding that 
(3.3) Vt > X t > . 

One may also use [3], p. 56, from which follows that (3.3) is equivalent to 

2aK 

— ^ 1 : 

a 1 

but, according to the above identifications, 

2an _ Z ~ A _ 2cf) _ 5 
er 2 a 2 a 2 

□ 

Our main result is the following : 
Theorem 3.4. Let us define the process 

z(t) = Jx t 

and the stopping time 

T = inf{t > 0\X t = 0}; 

as seen in Theorem 3.3, T = +oo a.s. for 5 > 2, and T < +oo a.s. for 6 < 2. 
Then there exists a Bernstein process y(t) for 

-I 
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and the potential 



V(t,q) = ^+Dq 2 . 



such that 

\Jte[0,T[z(t)=y(t) . 
In particular, for S > 2, z itself is a Bernstein process. 

Proof. One has (cf. Proposition 3.1 and its proof) 

dX t = ay/x~ t dw(t) + (acf> - \X t )dt 

(t)dw(t) + (a4> - Xz(t) 2 )dt . 
Taking now f{x) = y/x, we have 

1 // 1 _3 

Vx > / (x) = ^-j= and / (x) = --x * , 

therefore, for all t e]0,T[, f (X t ) = — ^— and f" (X t ) = -\z{t)-^. The application 
of Ito's formula now gives : 

dz(t) = d(f(X t )) 

= f\x t )dX t + \f"{X t ){dX t f 

1 (az(t)dw(t) + (a4> - \z{t) 2 )dt) - \z{t)~ s a 2 z(t) 2 dt 



2z(t) 

Let us now define r\ by 



= ^dw(t) + ^(4a4> ~ 4\z(t) 2 - a 2 )dt . 



\<j)t Xq 2 2<j> 1 
T](t,q) := e a a 2 q a 2 
XSt Xq 2 6-1 
= e~ ~ ~~tfq~^2~ . 

It is easy to check that 77 solves the equation 

2 dr, 9 4 d 2 n 

d-t=-YW + Vr] 



for 



in other words, 



C 

7 2 



V=^+Dq 2 



S := -6 2 ln(r)) 



16 4 

satisfies (HJB V ). Furthermore we have : 



n2 .XSt Xq 2 (5-1 
^ +^-o?{^)Hq) 
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B := 



whence 



drj 

V 
dq 

Ag a 2 (6-l) 
2 8q 

— (a 2 S - a 2 - 4\q 



B(t, z(t)) = g^y(4«0 - a 2 - 4Az(i) 2 



and z satisfies the stochastic differential equation associated with 77 : 

V* e]0, T[ dz{t) = 9dw{t) + B(t, z(t))dt 
(as in [5], §1, equation (B)) ; the result follows. □ 

Proposition 3.5. The isovector algebra T-Ly associated with V has dimension 6 if 
and only if 4> £ ^r}> i-e. S 6 {1,3} ; in i/ie opposite case, it has dimension 4. 



Proof. It is enough to apply Theorem 2.1, observing that the condition C = is 

a 3a 
4 ' 4 



equivalent to € {f , ¥}. □ 



In the context of Henon's already mentioned PhD thesis ([3], p. 55) we have 
4> = na, A = k, a = a 2 et /3 — 0, whence <j> — na and the condition C = is 
equivalent to 

r a 2 3a 2 , 

KflG{ T'~ } ' 

Let us analyze more closely the situation in which C = ; the general case will be 

commented upon in [5]. 
a 

1)0 = — , i.e. 6=1. 
Then y(t) is a solution of 

dy(t) = ~dw(t) - ~y(t)dt , 

i.e. y{t) is an Ornstein-Uhlenbeck process (it was already known that the Ornstein- 
Uhlcnbeck process was a Bernstein process for a quadratic potential). Therefore 
z(t) coincides, on the random interval [0, T[, with an Ornstein-Uhlenbeck process. 
Here 

Xt Xg 2 
T)(t,q) = e4 a 2 . 

From 

i/(t) = e-^(y Q + -J e^dwis)) 
_ xt / ~/a 2 (e At -l)^ 
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(w denoting another Brownian motion), it appears that y(t) follows a normal law 

l-e" 
4A 



with mean e ^ zq and variance " ^.f -. The density pt(q) of y(t) is therefore 



given by : 



2VX 2X(q~e-^z )\ 
Ptiq) = aV2,(l-e^) eXP{ - a'(l- e -*) } 



Whence 

Pt(q) 



Vt > 77* (*,£?) = 



-Xq 2 - Xq 2 e~ xt + 4Agz e"^ - 2Az 2 e~ 



1 



"A <" 



\| TTSinh(y) 



a 2 (l-e~ xt ) 



2 

and one may check that, as was to be expected, 77* satisfies the following equation 
(Cf } in 0) : 

2)0=^, ie. 5 = 3. 

In that case, according to Theorem 3.2, T = +00 whence y — z. Furthermore 



A „3a 2 i 
?y(f,g) = qea 

Let us define 

At 

s(t) = e'2±- ; 

z(t) 

then an easy computation, using Ito's formula in the same way as above, shows 
that 

ds(t) = -^e^s(t) 2 dw{t) ; 

in particular, s(t) is a martingale. 

Referring once more to Proposition 3.1 and its proof, we see that 

dX t = as/X~ t dw{t) + (^- - \X t )dt . 
Let us now assume Xq = and A 7^ ; then, according to Corollary 3.2, 

where Y is a BESQ 3 (squared Bessel process with parameter 3) such that Y(0) = 0. 
But, for each fixed t > 0,Y t has the same law as tY x , and Y\ = ||-E>i|| 2 is the square 
of the norm of a 3-dimensional Brownian motion ; the law of Y\ is therefore 



\=e %y/ul u > du 
2ir 
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Therefore the density pt(q) of the law of z(t) is given by : 

2\g 2 



Pt(q) = 



2ir a 3 (l — e~ A *)i 



and 



3Ai 



M 2 



Vi > rj*(t,q) 



Pt(q) 



16A5 



(1 



xtyi qe 4 a 2 t anh(f) 



Here, too, one may check directly that 77* satisfies equation (3.4) above. 
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